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Abstract — The progressive edge-growth algorithm is a well- 
known procedure to construct regular and irregular low-density 
parity-check codes. In this paper, we propose a modification of the 
original algorithm that improves the performance of these codes 
in the waterfall region when constructing codes complying with 
both, check and symbol node degree distributions. The proposed 
algorithm is thus interesting if a family of irregular codes with 
a complex check node degree distribution is used. 

Index Terms — Irregular low-density parity-check codes, pro- 
gressive edge-growth algorithm, waterfall region. 



I. Introduction 

Low-density parity-check (LDPC) codes were introduced 
by Gallager in the early 1960s [Ij. These codes are capacity 
achieving for many communication channels, and though 
forgotten for years, are nowadays ubiquitous as they have been 
found useful in many communications problems. 

A family of LDPC codes is usually defined by two gen- 
erating polynomials, X{x) and p{x). The coefficients of these 
polynomials define the distribution of incident edges to symbol 
and check nodes respectively. Richardson et al. showed that 
the asymptotic behaviour of a family of codes defined by both 
polynomials can be analysed by using the density evolution 
algorithm ||2l- Families of LDPC codes performing close to 
the channel capacity can be then designed by optimising both 
generating polynomials |3|, Q. 

An LDPC code is a linear code identified by an sparse 
parity-check matrix or its equivalent bipartite graph, also 
called Tanner graph. It is known that some iterative algorithms, 
such as beUef propagation based algorithms, provide optimum 
decoding over cycle-free Tanner graphs Q. However, any 
finite-length graph has necessarily cycles, and it has been 
shown that a large girth (length of the shortest cycle) improves 
the performance of LDPC codes using iterative decoding as 
it enforces a reasonable minimum distance [6|. Therefore, 
taking into account that a finite-length graph has cycles, it 
is important to make its girth as large as possible. 

The progressive-edge-growth (PEG) algorithm is an efficient 
method for constructing Tanner graphs with large girth ||6l, in 
most cases with better performance than randomly constructed 
codes. PEG algorithm's interest lies in its simplicity, and its 
flexibility when constructing codes from a complex symbol 
node degree distribution. However, note that a large girth does 
not automatically imply a a large minimum distance. The 
performance of these codes can be improved, for instance, in 
the error floor region. Significant research in PEG-based algo- 
rithms has been done to achieve this improved performance. 
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The rest of this paper has been organised as follows. In 
Section II, a modified PEG algorithm is proposed for the 
construction of LDPC codes following both X{x) and p{x) 
degree distributions for symbol and check nodes, respectively. 
In Section III, simulation results are shown. These results 
are compared with the original PEG algorithm and similar 
alternatives. In Section IV, the modified PEG algorithm and 
its applications in the construction of new codes are discussed. 

II. Modified PEG Algorithm 

A PEG-based algorithm consists of two basic procedures: a 
local graph expansion and a check node selection procedure. 
Both procedures are executed sequentially in order to construct 
a Tanner graph connecting symbol and check nodes in an 
edge-by-edge manner In the first procedure it is performed 
the expansion of the local graph from a symbol node, this 
expansion is used to detect and avoid short cycles when adding 
a new edge. The result is that check nodes that will produce a 
cycle are pruned, or if it is not possible to avoid a cycle, there 
only remains a set of candidate check nodes producing the 
largest cycle. The selection procedure is used to reduce this 
list of candidate nodes according to the current graph setting. 
In typical PEG algorithms, this procedure attempts to balance 
the degree of any check node selecting those candidates with 
the lowest check node degree. 

Recent research has been focused on the case when there are 
several candidate nodes after both procedures. At this point, 
it is possible to improve the performance of a PEG-based 
algorithm, for instance avoiding small stopping sets for the 
binary erasure channel (BEC) Q-lHl or trapping sets in the 
binary symmetric channel (BSC). New definitions have also 
been introduced, such as the extrinsic message degree (EMD) 
or the approximate cycle EMD (ACE) ifTOll . ifTTI . which are 
two common measures used to calculate the connectivity of 
symbol nodes. 

In the original PEG algorithm ||6l, a code is constructed 
according to a symbol node degree sequence. This sequence 
is previously calculated with the number of symbol nodes n 
and the edge degree distribution established by X{x). Note that 
the original proposal does not take into account the second 
polynomial, p{x), for the check node degree distribution. The 
algorithm proposed in this paper follows both degree distri- 
butions during the code construction procedures, changing the 
edge selection criterion, thus obtaining a better approach to 
codes with an irregular degree distribution. 

A. Free Check-node Degree (FCD) Criterion 

The edge selection procedure used in this proposal differs 
from the selection procedure proposed in the original PEG 
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Fig. 1 . Tanner graphs illustrating some characteristics of the algorithm. In 
(a) the zig-zag pattern used for 2-degree symbol nodes is shown, (b) and 
(c) show check nodes with different degrees, d{ci) and d{cj), and different 
number of incident edges (partially assigned), dJ'ici) = 1 and d^{cj) = 2. 
First check node, c^, is selected if a lowest check node degree criterion is 
used, while the sec 
(see Section Ill-A) . 



algorithm. The graph is analysed to avoid local short cycles, 
however check nodes are not chosen according to its number 
of assigned edges, d'^{ci), i.e. its current (or partial) degree. 

Instead of this, the check node with the highest differ- 
ence between its partial and final-defined degree is chosen, 
f{ci) — d{ci) ~d^{ci), i.e. the difference between the number 
of currently assigned edges and the total number of edges 
to be assigned. The lowest check node degree procedure is 
replaced by a highest free check node degree (FCD) procedure 
(see Fig. [TJ. The FCD concept, comes from the concept of 
"sockets" previously described in 121, lfT2l . 

We introduce the concept of compliance of a constructed 
code as the distance between the distribution of nodes (symbol 
or check nodes) in the code and the pre-established node 
degree distribution. Let pj be the pre-established probability 
distribution for a degree j check node, and p* the actual 
probability of a degree j check node in the constructed graph, 
we calculate the p-compliance of a code as: 



V=Y1 \P3-P*j\ 

J=2 



(1) 



else 

El. ■(r- {ci,Sj), where £'j. is the first edge 
incident to Sj and Ci is a check node such that 
it has the highest free check-node degree. 
end if 
else 

Expand a subgraph from symbol node Sj up to depth 
I under the current graph setting, such that A/"] = 

Mi+\ or ATI^. ^ but Arl+' = 0. 

E^. <— {ci,Sj), where E^. is the kth edge incident 

to Sj and Ci is a check node picked from the set 77^ 
having the highest free check-node degree. 
end if 

f{c^) = f{c^) - 1 

end for 
end for 

The same notation as in [6| is used, d{sj) is the Sj symbol 
node degree (The number of incident edges. It corresponds to 
the cardinality of the ensemble Eg after the code construc- 
tion), d{ci) is the Ci check node degree, /(q) is the number 
of edges that can be added to the check node Cj under the 
current graph setting, such that d*'{ci) = /(q) — d{ci). Eg- is 
the ensemble of edges connected to the symbol node Sj, E^. 
the edge added in the step k of the progressive construction, 
and Afg. is the ensemble of nodes reached after the graph 
expansion from the symbol node Sj up to depth /. 

A zig-zag construction for 2-degree symbol nodes (see 
Fig. [Til is forced by using a special criterion when adding 
the first edge to a symbol node. In this particular selection, a 
list of eligible check nodes is limited to those check nodes 
already connected under the current graph setting. Eg' = 
Eg^ U Eg^ U ■ • ■ U Eg._-^, i.e. to the fist of check nodes that 
have been chosen at least once from the first to j-th step. 
This construction is used to avoid cycles with 2-degree symbol 
nodes, thus obtaining a good performance in the error floor 
region as the results show (see Section Hill below). 

Relaxed edge selection: The proposed PEG algorithm can 
be modified to work with a relaxed edge selection. In this 
case, if there are not check nodes with free edges in the final 



where p*{x) — X]j=2 Pj^^ ' ^"'^ '^c^^^ — ma.x{d{cj)}. 



ensemble of candidate check nodes, 77 „ , check nodes with 
free edges are searched in the previous candidate ensemble, 
jVg . This procedure improves the p-compliance, r/, with the 
target check node degree distribution, p{x), at the expense of 
the current local cycle length (see Table |l|i. 



B. Modified Progressive Edge-Growth (PEG) Algorithm 

A modified PEG algorithm is described below. 
Require: d{si) < d{sj) Vi < j and /(c^) — d{ci) Vi 
for j = 1 to n do 
for k — 1 to d{sj) do 
if fc = 1 then 
i{ d{sj) = 2 then 

El. <— {ci,Sj), where El. is the first edge 
incident to Sj and Ci is a check node such that 
it has the lowest check-node degree under the 
current graph setting Eg.^ U Eg^ U • • • U Eg._-^ . 



III. Simulation Results 

Simulations results have been done for different LDPC code 
construction methods. All constructed codes have a codeword 
length n = 10^ and rate one half, R — 0.5. Performance has 
been measured under iterative decoding by using belief propa- 
gation. The maximum number of iterations for the decoder was 
set to 2000. Five different PEG-based construction methods 
are compared: (1) the original PEG algorithm as proposed 
in ||6J; (2) the modified PEG algorithm proposed by Richter 
in lfT2l : (3) a modification of (2) by inforcing the selection of 
a check node in the current graph when the first edge is added 
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TABLE I 
p-COMPLIANCE CALCULATED AS DEFINED IN EQ. (Q} 
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Fig. 2. Frame en'or rate (FER) as a function of the signal-to-noise ratio 
(SNR) for five construction methods, as explained in the text. Results have 
been obtained with the simulation of transmissions over the AWGN channel. 
Generating polynomials used in the code construction have been extracted 
from Table 11 (with a maximum symbol node degree of 50) in t2l . 
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Fig. 3. FER as a function of the crossover probability, e, over the BSC 
for five construction methods (see text). Codes have been constructed using 
the optimised generating polynomials from |4|. Since the p{x) distribution 
is more complex in this set of codes, it is possible to better appreciate the 
differences among the various algorithms used for constructing the codes. 



to a symbol node; (4) the modified PEG algorithm proposed 
here; and (5) a mixed version, the lowest check node degree 
criterion is used to connect the first edge to a symbol node 
(not only to 2-degree symbol nodes as proposed here) and the 
FCD criterion is used for the remaining edges. 

Figs. |2] and [3] show the performance of the codes over the 
additive white Gaussian noise (AWGN) channel and the BSC. 
The error floor is improved using the zig-zag construction for 
2-degree symbol nodes. On the other hand, within a given 
graph setting, when the first check node connected to a symbol 
node rule is used, there is no relevant improvement, as this 
allows for different structures to the zig-zag. 
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show that this method behaves very well in the waterfall region 
while also maintaining a low error floor. Good performance 
in both regions is relevant as it allows to use the same code 
independently of the working point, e.g. this construction 
might find application in rate compatible solutions. 
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IV. Conclusions 

In this paper, we proposed an improved PEG algorithm. 
The constructed codes comply with both an irregular check 
and an irregular symbol degree distribution. Simulation results 
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